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Computational  Study  of  Chalcopyrite 
Semiconductors 

and  Their  Non-Linear  Optical  Properties 


Walter  R.  L.  Lambrecht 

1  Publications 

a.  Publications  published  in  peer  reviewed  journals 


1.  “Electronic  band  structure  of  ordered  vacancy  defect  chalcopyrite  com¬ 
pounds  with  formula  Xiaoshu  Jiang  and  Walter  R.  Lam¬ 

brecht,  Phys.  Rev.  B  69.  035201  (2004) 

2.  “Noncritically  phase-matched  second-harmonic-generation  chalcopy- 
rites  based  on  CdSiAs2  and  CdSiP2,"  Walter  R.  L.  Lambrecht  and 
Xiaoshu  Jiang,  Phys.  Rev.  B  70.  045204  (2004) 

3.  “Theoretical  study  of  cation- related  point  defects  in  ZnGeP2,”  Xiaoshu 
Jiang,  M.  S.  Miao,  and  Walter  R.  Lambrecht,  Phys.  Rev.  B  71,  205212 
(2005). 

4.  “Structure  and  phonons  of  ZnGeN2,"  Walter  R.  Lambrecht,  Erik  All¬ 
dredge,  and  Kwiseon  Kim  Phys.  Rev.  B  72,  155202  (2005) 

5.  “Theoretical  study  of  the  phosphorus  vacancy  in  ZnGeP2,”  Xiaoshu 
Jiang,  M.  S.  Miao,  and  Walter  R.  Lambrecht  Phys.  Rev.  B  73.  103203 
(2006) 

G.  “Large  band-gap  bowing  in  C1i1_2.Agj.GaS2  chalcopyrite  semiconduc¬ 
tors  and  its  effect  on  optical  parameters,"  Chandrima  Mitra  and  Wal¬ 
ter  R.  Lambrecht,  Phys.  Rev.  B  76,  035207  (2007) 

7.  “Universal  Transition  State  for  High-Pressure  Zinc  Blende  to  Rocksait 
Phase  Transitions,”  M.  S.  Miao  and  Walter  R,  Lambrecht,  Phys,  Rev. 
Lett.  94,  225501  (2005) 
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b-  Manuscripts  submitted  to  peer  reviewed  journals  but  not  yet,  published 


1.  “First-principles  calculation  of  the  zone  center  phonons  in  ZuSiNo: 
comparison  with  infrared  data,'1  Tula  R.  Paudel  and  Walter  R.  L. 
Lambrecht,  Phys.  Rev.  R  (2007),  accepted,  in  press. 

2,  “Electronic  and  crystal  structure  of  Cua-aS,”  Pavel  Lukashev,  Walter 
R,  L.  Lambrecht,  Takao  Kotani  and  Mark  van  Schilfgaarde,  submitted 
to  Phys.  Rev.  B  (July  2007) 

3*  “Baud  gap  bowing  in  AgGa(SeI_^Te^)2  system  and  its  effect  on  the 
second  order  response  coefficient  and  the  refractive  index,”  Chandrima 
Mitra  and  Walter  R*  L.  Lambrecht,  submitted  to  Phys.  Rev  B  (July 
2007) 

4.  “First-principles  study  of  the  oxygen  vacancy  in  ZuO:  self-consistent 
gap-corrected  approach,”  Tula  R Paudel  and  Walter  R.  L.  Lambrecht, 
submitted  to  Phys.  Rev  B  (August  2007), 

c.  Publications  in  conference  proceedings 

1.  “Electronic  Structure  of  Native  Point  Defects  in  ZnGeP^,”  Xiaoshu 
Jiang*  Maosheng  Miao  and  Walter  R*  L.  Lambrecht,  in  Progress  in 
Compound  Semiconductor  Materials,  Ill-Electronic  and  Optoelectronic 
Applications ,  ed.  Daniel  J.  Friedman.  Omar  Manasreh*  Irina  A*  Buyanova, 
Amieli  Munkholm,  F.  Danie  Auret,  Mat,  Res,  Soc,  Symp.  Proc,  VoL 
799*  Z5.3.1  (2004)* 

d.  Presentations  at  conferences  not  published  in  proceedings 

L  "First-principles  study  of  native  point  defects  in  ZnGePj,”  Xiaoshu 
Jiang  and  Walter  R.  L.  Lambrecht,  Fall  Meeting  of  the  Ohio  Sec¬ 
tion  October  17-18,  2003,  Case  Western  Reserve  University  Cleveland, 
Ohio*  Abstract  P.005 

http://flux.aps.org/meetings/YR03/OSF03/baps/abs/Sl(K)0l,html 

2,  “Point  defects  and  defect  complexes  in  ZnGeP2<v  Xiaoshu  Jiang,  Maosheng 
Miao  and  Walter  R*  L.  Lambrecht,  APS  March  meeting  2004  (Mon¬ 
treal)*  abstract  K  1.037  Bulb  Am.  Phys.  Soe*  49  (2004); 
http;//flux*aps*org/meetings/YR04/MAR04/baps/abs/S3610037Jitml 
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3*  “First-principles  Calculations  Based  Desing  of  Chalcopyrite  Semicon¬ 
ductors  for  Nonlinear  Optical  frequency  Conversion,”  Walter  R.  L. 
Lambrecht,  Enabling  Photonics  Technologies,  Great  Lakes  Photonics 
Symposium  June  7-11,  2004,  Cleveland 

4,  “Electronics  structure  of  simple  and  complex  point  defects  in  ZnGeP2,1T 
X.  Jiang,  M.  S,  Miao  and  W.  R.  L.  Lambrecht,  Research  Showcase 
2004,  at  Case  Western  Reserve  University,  April  2,  2004, 

5-  “Does  the  zinc  vacancy  in  ZnGePs  exhibit  a  Jahn- Teller  distortion?” 
Walter  R.  L.  Lambrecht  and  Xiaoshu  Jiang,  APS  March  Meeting,  Los 
Angelas,  March  2005,  Bull.  Am.  Phvs.  Soc.  50,  abstract  P  18.00008: 
http://meetings.aps.org/Meetmg/MARjD5/Gvent/23534 

6,  “Interactions  between  native  point  defects  in  ZnGeP2,”  Xiaoshu  Jiang, 

M.  S.  Miao,  and  Walter  R.  L.  Lambrecht,  APS  March  Meeting,  Los 
Angeles,  March  2005,  Bull,  Am.  Phys.  Soc.  50  Abstract  PI  8.00005; 
http://meetings,aps.org/Meeting/MAR(}5/Event/23530 

7,  “Electronic  bandstructure,  crystal  structure  and  phonons  of  ZnSiN^,1' 
Tula  R.  Paudel  and  W.  R,  L.  Lambrecht,  APS  March  Meeting,  Balti¬ 
more  2006,  abstract  R4 6,00009, 

http://meetings.aps.org/Meeting/MAR06/Event/44174 

8,  “Structural  and  Electronic  Properties  of  Copper  Sulfides,”  P.  Lukashev 
and  W.  R.  L,  Lambrecht,  18th  Annual  Workshop  on  Developments  in 
Electronic-Structure  Methods,  Columbus,  OIL  22-25  June  2006 

0.  “Structure,  electronic  structure  and  phonons  of  Zn-IV-N2  (IV=Sl,Ge,Sn),” 
T.  R.  Paudel  and  W.  R.  L.  Lambrecht,  poster  presented  by  T.  Rt 
Paudel  at  the  summer  school  on  Ab-Initio  Simulation  of  Crystalline 
Systems  2006,  Spokane,  WA  September  17-22,  2006. 

10.  “Point  defects  and  nonlinear  optics  modeling  in  chalcopyrite  semicon¬ 
ductors/7  Walter  R.  L.  Lambrecht,  seminar  at  National  Renewable 
Energy  Laboratory,  June  17  2005 

11.  “Point  defects  in  ZnGeP2,17  Walter  R.  L.  Lambrecht,  Condensed  Mat¬ 
ter  Seminar,  February  27,  2006,  Department  of  Physics,  Case  Western 
R  eser  ve  Uni  vers  i  ty 

12.  “Large  band  gap  bowing  in  CuAgGaS2  chalcopyrite  semiconductors 
and  its  effect,  on  optical  parameters,”  C.  Mitra  and  W.  R.  L,  Lam¬ 
brecht,  APS  March  Meeting,  Denver  2007,  Bull.  Am.  Phys,  Soc,  52, 
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abstract  no.  D40.0GQQ8  (2007); 
http://ineetings.aps.org/Meeting/MAIl07/Event/57971 

13.  “Crystal  and  Electronic  Structure  of  Copper  Sulfides,”  P>  Lukas  he  v 
and  W.  R.  L.  Lambrecht,  APS  March  Meeting,  Denver  2007,  Bull. 
Am.  Phys.  Soc.  52,  abstract  no,H4L00013  (2007); 
http://meetings.aps.org/Meeting/MAR07/Event/58515 

14.  “Electronic  structure  of  the  oxygen  vacancy  in  ZnO  including  a  self- 
consistent  band-gap  corrected  approach,”  Walter  R.  L  Lambrecht  and 
T,  R.  Paudel,  to  be  presented  at  the  MRS  Fall  Meeting  2007  in  Sym¬ 
posium  L:  ZnO  and  Related  Materials,  abstract  L4.2. 

15.  “Phonons  in  Zn-IV-N2  semiconductors,”  T.  R.  Paudel  and  W.  R.  L. 
Lambrecht,  to  be  presented  at  the  MRS  Fall  Meeting  2007  in  Sympo¬ 
sium  Q:  Nitrides  and  Related  Bulk  Materials,  abstract  QL9. 

1G.  “Raman  Spectroscopy  of  Single  Crystal  ZnGeMa,”  Timothy  Peshek, 
Kathleen  Kash,  John  C.  Angus,  Tula  R.  Paudel  and  Walter  R.  L. 
Lambrecht,  to  be  presented  at  the  MRS  Fall  Meeting  2007  in  Sympo¬ 
sium  Q:  Nitrides  and  Related  Bulk  Materials,  abstract.  Ql.l 

e.  Manuscripts  in  preparation  to  be  submitted  to  peer  reviewed  journals 

1.  “Jalm-Teller  distortion  of  the  zinc  vacancy  in  ZnGePs,”  Xiaoshu  Jiang 
and  Walter  R.  L,  Lambrecht,  in  preparation 

2.  “Native  point  defects  in  CdGeAss,”  Tula  R.  Paudel  and  Walter  R.  L. 
Lambrecht,  in  preparation 

L  Theses 


1.  “First  principles  study  of  electronic  structures  of  defects  in  ZnGeP2 
and  defect  chalco pyrites,”  Xiaoshu  Jiang,  Ph,  D,  Thesis,  Case  Western 
Reserve  University  (2005),  available  from  Kelvin  Smith  Library. 

2.  “Crystal  and  electronic  structure  of  copper  sulfides”,  Pavel  Lukashev. 
Ph,  D,  Thesis,  Case  Western  Reserve  University  (2006),  available  from 
Kelvin  Smith  Library  and  in  electronic  form: 
http://www.ohiolink.edu/etd/view.cgi7acc 

Summary 

a.  Number  of  published  papers  in  peer  reviewed  journals:  7 
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b*  Number  of  submitted  (not  yet  published)  papers  to  peer  reviewed  jour¬ 
nals:  4 

c.  Number  of  conference  proceedings  published:  1 

d*  Number  of  presentations  at  conferences,  seminars,  not  published:  16 

e.  Number  of  manuscripts  in  preparation  (not  yet  submitted):  2 

f.  Number  of  theses  completed:  2 


2  Summary  data  on  human  resources 

The  following  scientists  were  involved  with  the  project: 

•  Walter  R,  L*  Lambreeht,  Principal  Investigator,  Professor  of  Physics 

•  Maosheng  Miao,  postdoctoral  Sr.  Research  Associate  (part  time,  2002- 
2006) 

•  Xiaoshu  Jiang,  graduate  student  (full  time,  completed  PIlD,  2002- 

2005) 

•  Pavel  Lukashev,  graduate  student  (part  time,  completed  Pli.D*  2005- 

2006) 

•  Tula  Ram  Paudel,  graduate  student  (part  time,  PIlD.  in  progress, 
2005-2007) 

»  Chandrima  Mitra,  graduate  student  (part  time,  Ph.D*  in  progress, 
2005-2007) 

a.  Number  of  Scientists  supported  by  this  agreement:  6 

b.  Number  of  Graduate  Students  supported  by  this  agreement:  4 

c.  Number  of  Post.  Doctoral  Associates  supported  by  this  agreement:  I 

d.  Number  of  Faculty  supported  by  this  agreement:  1 

e.  Number  of  Other  Staff  supported  by  this  agreement;  0 

f.  Number  of  Undergrads  supported  by  this  agreement:  0 

g-  Number  of  PhD(s)  awarded  as  a  result  of  this  agreement:  2 

h.  Number  of  Bachelors  Degrees  awarded  as  a  result  of  this  agreement:  0 

i.  Number  of  Master  Degrees  awarded  as  a  result  of  this  agreement:  0 

j.  Number  of  FTE  graduate  students  supported  per  year:  1 

k.  Number  of  FTE  postdoctoral  associates  supported  per  year:  I 
L  Faculty  summer  salary  support  in  months  per  year:  1 
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Report  of  Inventions 


a.  Number  of  inventions  or  discoveries  disclosed  to  University  Office  of  Spon¬ 
sored  Research:  1 

Title:  Materials  that  allow  for  lion-critical  phasematching  in  second  order 
nonlinear  optical  frequency  conversion 
Date  of  disclosure:  September  14,  2007. 

b.  Number  of  Patents  Sumbitted  as  a  result  of  this  agreement:  0 

c.  Number  of  Patents  Awarded  as  a  result  of  this  agreement:  0 


4  Narrative  Description  of  Scientific  Accomplish¬ 
ments 

The  main  goals  of  the  project  were  to  study:  (1)  the  physics  of  point  defects 
in  ehalcopyrite  semiconductors  that  have  currently  applications  in  nonlinear 
optics,  such  as  ZnGePs  and  CdGeAs^;  (2)  the  properties  of  ternary  alloys 
among  different  ehalcopyrite®  to  evaluate  the  dependence  of  optical  proper¬ 
ties  on  the  composition  of  the  materials;  (3)  understand  the  trends  in  the 
family  of  ehalcopyrite  semicondutors  and  related  materials  with  the  aim  to 
identify  potential  materials  with  favorable  combinations  of  properties  for 
non-linear  optical  applicsations. 

The  methodology  followed  in  the  work  is  based  generally  on  density 
functional  theory  first-principles  calculations.  We  use  the  full-potential  lin¬ 
earized  muffin-tin  orbital  method  to  solve  the  Kohn-Sham  equations  in  the 
local  density  approximation  (LDA).  Defects  are  modeled  with  supercells  (Le. 
peridodic  boundary  conditions).  Alloy  problems,  in  particular  the  disorder 
aspect,  were  tackled  with  the  use  of  so-called  special  quasirandom  struc¬ 
tures  (SQS).  Several  aspects  of  the  methodology  were  refined  during  the 
course  of  the  project.  For  example,  early  on,  we  had  to  include  the  use  of 
a  uniform  background  compensating  charge  density  in  the  codes  in  order  to 
allow  us  to  deal  with  charged  point  defects.  The  issue  of  size  dependence 
in  such  calculations,  and  so-called  Makov-  Payne  monopole  and  quadrupole 
correctors  were  studied.  We  also  demonstrated  the  use  of  shift  potential  to 
address  the  LDA  band  gap  problem  in  point  defect  studies  and  ultimately 
showed  how  this  problem  can  be  more  consistently  be  addressed  by  a  novel 
implementation  of  the  LDA+U  approach.  Thus,  several  improvements  in 
methodology  resulted  from  this  work. 

The  following  tasks  were  accomplished  and  led  to  the  publications  iden- 
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titled  below, 

•  A  study  was  carried  out  of  the  cation  antisite  defects  Zuqg,  Gezn, 
cation  vacncies  V^nt  Vqg  and  anion  vacancies  Vp  in  ZnGePs  and  some 
complexes  among  these  defects  and  resulted  in  3  publications  [1,  2,  3] 
and  a  major  portion  of  the  Ph,  D.  thesis  of  Xiaoshu  Jiang. 

•  A  further  in-depth  study  was  carried  of  the  possibility  of  a  Jahn- 
Teller  distortion  of  the  Vzn  defect  .  While  a  possible  mechanism  which 
could  explain  the  associated  experimental  {electron  paramagnetic  res¬ 
onance  and  electron  nuclear  ddUuble  resonance)  EPR-ENDOR  data 
was  identified, [4]  some  technical  problems  remain  associated  with  this 
problem.  The  currently  available  level  of  theory,  which  is  based  on 
the  local  density  approximation  (LDA),  does  not  find  the  proposed 
Jalm- Teller  distortion  to  be  energetically  favorable.  It  was  found  that 
in  part  this  could  be  related  to  the  formation  of  Ge  dangling  bonds 
next  to  the  vacancy  as  a  result  of  using  too  small  a  simulation  cell  but 
also  could  be  related  to  the  incomplete  cancellation  of  self-interaction 
effects  in  LDA.  A  calculation  at.  the  Hartree-Fock  level  was  pursued 
using  duster  models  and  the  Gaussian  program  but  did  not  lead  to  a 
satisfactory  explanation  of  the  experimental  results,  A  description  of 
this  work  is  included  in  the  Ph.  D,  Thesis  of  Xiaoshu  Jiang, 

•  A  similar  study  was  carried  out  of  the  point  defects  in  CdGeAss-  Be¬ 
cause  CdGeAs2  has  a  negative  band  gap  in  the?  local  density  approx¬ 
imation,  it  was  essential  to  develop  a  band  gap  correction  approach. 
This  was  done  using  the  so-called  LSDA+U  method,  which  was  de¬ 
veloped  under  another  project  (supported  by  ONR  and  NSF)  during 
the  period  of  this  grant.  The  results  were  only  obtained  in  the  last 
months  of  the  extension  period  of  this  project  and  remain  to  be  fully 
analyzed  and  prepared  for  publication, 

•  As  part  of  our  tests  of  the  LSDA-bU  approach  for  dealing  with  band 
gap  corrections,  we  applied  the  method  to  a  recent  problem  of  in¬ 
terest  in  wide-band-gap  semiconductors,  namely  the  oxgygen  vacancy 
in  ZnO,  For  this  system,  a  controversy  exists  in  the  literature  as  to 
whether  the  defect  levels  lie  in  the  uppor  or  lower  half  of  the  band 
gap,  and  this  discrepancy  results  purely  from  different  ways  of  apply¬ 
ing  band  gap  corrections.  Our  new  approach  sheds  new  light  on  this 
problem,  A  paper  on  this  work  was  submitted  for  publication  [5]  and 
the  work  was  also  submitted  for  presentation  at  the  upcoming  MRS 
Fall  meeting. 
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*  A  study  was  carried  out  of  the  electronic  band  structures  in  the  family 
of  II-IIWLi  defect  chalcopyrites,  These  are  related  to  l-III-Vl^  chal- 
copyrites  by  a  doubling  of  the  unit  cell  and  replacing  the  group-I  atoms 
by  a  group  If  atom  and  a  vacancy.  These  materials  have  received  only 
little  attention  in  the  literature  and  our  work  [G]  presents  the  first  sys¬ 
tematic  approach  to  understanding  the  trends  in  their  band  gaps  and 
other  properties. 

*  In  terms  of  optimizing  materials  for  non-linear  optical  calculations, 
we  had  found  in  a  previous  funding  period  supported  by  AFOSR,  that 
AgGaTez  has  a  twice  higher  \ ^  than  AgGaSe-2,  which  is  currently  one 
of  the  main  optical  materials  used  for  frequency  conversion  because  of 
its  very  wide  window  of  transparency.  Work  at  WPAFB  by  Mel  Ohmer 
and  coworkers,  had  also  found  that  AgGaT©2  has  a  positive  instead 
of  negative  birefringence.  This  lead  to  the  idea  (which  was  in  fact 
patented  by  researchers  at  WPAFB)  that  alloys  of  AgGa(Sei_3:Te3T)2 
could  be  designed  to  have  the  optimal  index  of  refraction  to  allow  for 
so-called  non-critical  phasematching.  However,  this  all  rests  on  the 
assumption  that  indices  of  refraction  vary  linearly  with  concentration 
in  the  alloy*  However,  usually  semiconductor  band  gaps  show  some  de¬ 
gree  of  bowing  instead  of  linear  behavior*  In  turn  this  can  he  expected 
to  lead  to  a  nonlinear  behavior  of  the  indices  of  refraction.  We  thus 
studied  the  fundamental  question  of  how  the  index  of  refraction  and 
the  vary  as  function  of  concentration  in  quaternary  alloy  systems. 
We  carried  out  this  study  for  AgGatSei^Te^  in  the  last  months  of 
the  project  and  submitted  the  results  for  publication* [7]  It  was  found 
that  in  this  system  the  bowing  of  the  index  of  refraction  is  only  mi¬ 
nor,  so  that  previous  results  based  on  linear  behavior  should  provide 
a  reasonable  estimate* 

*  Closely  related  to  the  previous  topic,  we  also  studied  the  band  gap 
bowing  and  index  of  refraction  and  as  function  of  concentration 
x  in  Cu^Agi-tf  GaS-2  alloys*  For  this  quaternary  alloy  family,  exper¬ 
imental  reports  existed  of  a  large  band  gap  bowing.  We  found  that 
this  leads  indeed  to  large  bowing  of  the  index  of  refraction  and  of 

An  interesting  feature  of  this  alloy  system  is  that  while  Cu-based 
compounds  have  smaller  gaps  than  the  corresponding  Ag- based  com¬ 
pounds,  they  nevertheless  have  a  smaller  in  defiance  of  the  rule  of 
thumb  that  smaller  gaps  lead  to  higher  that  is  generally  observed. 
In  fact,  an  unexpected  result  of  this  work  was  that  the  °f  a  50  % 
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mixed  quaternary  Cuo^Ago^GaSg  is  predicted  by  our  calculations  to 
have  a  higher  than  either  of  the  end  compounds.  This  work  was 
published  in  [8], 

•  When  comparing  the  two  families  of  chalcopyrite  materials,  the  II- IV- 
V2  family  and  the  I-III-VI2  family,  one  generally  finds  higher  and 
also  stronger  bonds  in  the  former  than  in  the  latter.  Stronger  bonds  are 
among  other  expected  to  lead  to  higher  laser  damage  thresholds,  and 
are  thus  important  for  improving  the  tolerable  input  and  output  power 
levels  of  frequency  conversion  optical  devices.  On  the  other  hand,  the 
birefringence  in  most  II-TV-V^  materials  is  positive,  thus  precluding 
so-called  non-critical  phase  matching.  This  led  us  to  the  question:  do 
there  exist  II-IV-V2  chalcopyrites  with  negative  birefringence,  and  if 
so  can  they  be  alloyed  to  achieve  non-critical  phasematching?  As  al¬ 
ready  outlined  in  the  proposal  for  this  project,  we  found  a  correlation 
between  birefringence  and  c/a  ratio  of  the  crystal  structure  and  on  this 
basis  anticipated  that  CdSiAs2  and  Cd  SiP-2  would  have  negative  bire¬ 
fringence.  We  thus  embarked  on  a  computational  study  of  these  ma¬ 
terials.  We  found  indeed  that  both  have  negative  birefringence  by  our 
calculations,  and  in  fact,  later  found  confirmation  of  this  in  some  little 
known  works  in  the  Russian  literature,  which  was  not  yet  included  in 
the  reference  compilations  on  which  our  original  proposal  was  based. 
However,  when  studying  in  more  detail  whether  non-critical  phase- 
matching  is  possible,  for  some  desirable  target  frequency  conversion 
processes  into  the  mid  infrared  range,  we  found  that  the  birefringence 
in  CdSiAss  is  too  small  to  overcome  the  index  of  refraction  disper¬ 
sion  between  pump  and  target  beam  wavelengths.  However,  our  work 
also  suggests  that  by  compressing  the  c/a  ratio,  one  might  make  the 
birefringence  more  negative.  The  details  of  the  phasematching  angle 
curves  as  function  of  strain  were  worked  out  and  led  to  a  new  pro¬ 
posal  for  turning  CdSiAsg  into  a  useful  nonlinear  optical  material  bv 
applying  external  uniaxial  or  biaxial  stress.  In  fact,  this  concept  was 
ultimately  submitted  as  an  invention  disclosure.  CdSiPj  was  found  to 
have  stronger  negative  birefringence,  but  on  the  other  hand  has  a  less 
strong  So,  as  usual  there  is  a  trade-off.  One  could  tune  birefrin¬ 
gence  by  making  alloys  of  the  type  CdSifPi^As*^  This  work  was 
published  in  [9].  The  practical  implementation  of  these  ideas  will  how¬ 
ever  take  substantial  further  development  of  crystal  growth  techniques 
and  is  thus  not  close  to  commercial  applications. 
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•  Similar  to  the  chalcopyrite  family  which  is  based  on  an  ordered  ar¬ 
rangement  of  II  and  IV  cations  on  the  cation  sublattice  of  zincblende, 
there  exists  a  family  based  on  the  wurtzite  structure,  which  occurs 
naturally  for  111-nitride  materials.  These  are  potentially  an  interest¬ 
ing  family  of  optoelectronic  materials  competing  with  the  Ill-nitrides. 
We  had  in  the  past  already  studied  some  of  these  nitrides  but  as  an 
exploratory  subproject  of  this  grant,  we  started  a  study  of  the  lattice 
vibrations  in  these  materials.  We  did  so  mainly  because  we  were  aware 
of  the  availability  of  some  Raman  and  infrared  data  in  the  literature. 
During  a  sabbatical  stay  in  2005  at  NR  EL,  we  became  familiar  with 
computational  methods  for  calculating  phonons  in  crystals,  based  cm 
the  linear  response  method  and  the  A  BINIT  code.  This  led  to  a  first, 
publication  on  the  phonons  in  ZiiGeN^II]  In  this  work  however ,  we 
could  only  compare  with  Raman  spectra  on  polyerystalline  powders. 
Subsequently,  we  applied  the  same  approach  to  ZnSiN2  and  used  it 
to  simulate  infrared  reflectivity  spectra  on  single  crystals.  This  led  to 
a  second  publication,  [12]  We  are  now  trying  to  find  funding  for  con¬ 
tinuing  this  project  in  collaboration  with  an  experimentalist  in  our 
department  who  is  working  on  the  crystal  growth  of  this  family  of  ma¬ 
terials.  A  joint,  experimental/ computational  study  w as  carried  out.  of 
the  Raman  spectra  of  single  crystal  ZnGeNa  needle  type  crystals  and 
will  be  presented  at  the  MRS  Fall  Meeting  in  November  2007,  (See 
presentations  listed  above). 

*  Chalcopyrite  semiconductors  of  the  family  Cu(In,Ga)Se2  (CI(G)S)  are 
an  attractive  material  for  photovoltaic  (PV)  solar  energy  conversion 
in  thin  film  materials.  In  part,  this  is  owing  to  a  relatively  greater  tol¬ 
erance  to  lattice  defects  in  these  materials.  For  example,  high  conduc- 
tiviy  and  reduced  recombination  along  grain  boundaries  is  observed  in 
some  of  these  films.  This  in  turn  was  tentatively  related  to  the  nature 
of  the  valence  band  maximum  which  is  Cu*d  like.  This  leads  to  a  fa¬ 
vorable  band-offset  for  carrier  separation  at  interfaces  or  surfaces  with 
reduced  Cu  coordination.  However,  one  of  the  things  that  is  slowing 
down  the  development  of  these  PV  materials  is  the  relative  scarcity 
of  In,  It  thus  seemed  useful  to  look  for  other  Cu-based  compounds 
without  In.  We  found  through  literature  search  that  CujS  is  one  of 
the  earliest  used  solar  cell  materials  but  was  eventually  abandoned  be¬ 
cause  of  long-term  stability  problems.  We  thought,  this  question  was 
worth  re-investigating.  Encouraged  by  the  prospective  of  a  possible 
collaboration  with  a  chemist  in  our  university  who  was  able  to  syn- 
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thesize  CihS  in  nanoparticle  form,  wc  thus  started  a  study  of  Cu2S.  It 
turned  out  that  in  spite  of  its  simple  chemical  composition,  its  crystal 
structure  is  quite  complex  and  theoretical  studies  of  its  hand  structure 
had  never  been  reported  prior  to  ovir  study.  A  novel  methodology  had 
to  be  developed  to  construct  suitable  models  for  simulating  the  ran¬ 
dom  distribution  of  Cu  atoms  in  this  material.  Furthermore,  it  proved 
important  to  include  corrections  beyond  LDA.  To  this  end  we  collabo¬ 
rated  with  Mark  van  Schiifgaarde  at  ASU  who  provided  us  with  some 
calculations  using  his  recently  developed  state-of-the  art  quasiparticle 
self-consistent  GW  method  For  one  of  the  simpler  model  structures. 
We  were  then  able  to  use  this  information  to  apply  a  correction  to 
the  more  complex  models.  Our  work  on  this  subject  ultimately  led  to 
the  doctoral  thesis  of  one  Ph,D  student  and  one  paper  submitted  to 
Phys.  Rev*  B.[1G]  We  also  used  it  as  a  starting  point  for  proposals 
to  NSF  and  DOE  in  which  the  PV  applications  of  these  materials  in 
nanoparticle  form  would  be  pursued. 

Summary  and  outlook 

All  in  all,  most  of  the  initially  set  goals  of  this  project  were  successfully  ac¬ 
complished.  There  remains  work  to  be  done  on  point  defects  in  the  AgGaTej 
and  AgGaSes  systems  and  we  would  be  happy  to  pursue  this  with  future 
funding  from  AFOSR  if  there  is  a  continuing  interest  in  this  topic.  There 
also  remains  work  to  be  done  to  further  refine  the  computational  techniques 
for  dealing  with  some  of  the  more  subtle  aspects  of  defect  physics,  in  par¬ 
ticular  relating  to  the  partially  unresolved  problem  of  the  structural  relax¬ 
ation  around  a  Vzn  in  ZnGePg  and  the  associated  details  of  the  hyperfme 
structure  measured  in  magnetic  resonance  experiments.  Similar  problems, 
in  fact,  plague  other  known  defect  systems,  although  they  have  not  been 
widely  recognized  or  acknowledged.  The  project  led  to  several  spin-offs  into 
new  directions:  (1)  point  defects  in  ZuO,  (2)  the  fundamental  properties  of 
ternary  nitrides  of  the  type  II-IV-Ns,  a  new  family  of  wide-band-gap  semi¬ 
conductors,  (3)  the  study  of  Cu-sulfides  and  their  potential  in  photovoltaic 
energy  conversion* 


5  Interactions 

During  the  course,  of  the  project,  we  engaged  in  several  collaborations.  Our 
work  on  the  ZnGePj  point  defects  was  closely  related  to  the  EPR-ENDOR 
studies  by  the  group  of  Larry  Halliburton  and  Nancy  Giles  at  West-Virginia 
Univeristy,  Morgantown,  WV.  Our  work  on  ZnGeN2  and  related  materials, 
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was  initially  carried  out  with  researchers  at  NREL  (kwiseon  Kim)  and  is 
now  continued  with  a  collaborator  (Kathleen  Kasli)  at  CWRU.  Our  work  on 
C112S  was  inspired  in  part  by  work  of  a  colleague  in  the  Dept,  of  Chemistry 
(Clemens  Burda)  at  CWRU.  Our  work  on  (Cu,Ag)GaS2  was  inspired  by  and 
relevant  to  experimental  work  by  Prof.  Peter  Yu  at  University  of  California, 
Berkeley  and  his  collaborators  in  South  Korea. 

6  Summary  of  timeline  and  financial  report 

The  project  stretched  over  a  somewhat  unusually  long  period,  from  2002  t  ill 
2007.  It  was  initially  a  3  year  continuation  grant  for  a  total  of  $382,125, 
following  on  an  earlier  related  AFOSR  project  started  in  1999  and  extended 
with  a  no-cost  extension  till  Oct.  2002.  The  present  project  was  originally 
scheduled  to  end  on  October  14.,  2005,  but  received  a  no-cost  extension  till 
June  2006.  In  May  2006,  a  final  extension  with  additional  support  to  the 
amount  of  $30,000  was  requested  and  awarded  on  August  10,  2006  bringing 
the  total  to  $412,125.  The  new  extension  was  finally  no-cost  extended  till 
April  14,  2007.  The  present  final  project  report  was  delayed  by  the  desire 
to  first  finalize  some  of  the  papers,  so  as  to  enable  us  to  prepare  a  more 
accurate  final  report  representing  the  final  outcomes  of  the  project.  The 
writing  of  5  of  the  reported  papers  was  completed  during  the  summer  of 
2007,  after  completion  of  the  nffieal  grant  period. 
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